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Abstract 

We report on a confined quark model calculation of the transition form factors 
for the radiative decays of 77(77') — ► and 77(77') — ► jtt + tt^. The 77/77' mixing 

angle Op is determined from the set of data on the electromagnetic 77(77') decays. It 
is found that Op ~ —16.5°. 

The analysis of the dipionic mass distribution in the decay 77' — ► 77r + 7r~ confirms 
the existence of a non-resonant contribution which is the box anomaly. 
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1 Introduction 



The interest of studying the properties of the neutral pseudoscalar particles 7r°, r/ and 
T]' is not weaken from the begining of creation of unitary SU(3) symmetry. One of the 
intrigued problem was the so-called "U(l)-problem" which was originally encountered by 
Glashow 0. He considered extending SU(3)x SU(3) to U(3)x U(3) and found that if 
the extra U(l) axial symmetry was realized spontaneously, there would have to be an 
extra light pseudoscalar meson with mass compatible with pion (140 MeV). But the only 
flavor singlet meson is rf with a mass being equal to 958 MeV! After the advent of QCD, 
one could suggest that the anomaly associated with axial U(l) transformations might 
be used to provide a solution of U(l)-problem because of the divergence of axial current 
does not vanish in the chiral limit. However, the matter was found to be not so simple 
(see, for review ||). Witten || found that the anomaly could be turned off as a number 
of quark colors goes to infinity (N c — > oo) and then the 7/-meson would be considered 
as the Nambu-Goldstone boson. With finite N c he suggested that the mass squared of 
the f]' was proportional to 1/N C . Following this way, Veneziano ||J] proposed a more 
explicit mechanism which was consistent with anomalous Ward identities and realized 
it by using the effective chiral Lagrangian incorporated explicitly {7^(1) anomaly. The 
phenomenological implications of those Lagrangians to the 7] — ^'-mixing allowed one to 
get the mass formulas for the masses of r\ and v( and the mixing angle which reproduce 
the experimental values with quite resonable accuracy. 

The three-flavor Nambu-Jona-Lasinio (NJL) model with the instanton-induced six- 
quark effective interaction which incorporates effects of {7a (1) anomaly has been developed 
to study the spectrum of low-lying mesons, the r\ — //-mixing, etc. ||. 

The lattice calculations of the rf mass have been performed in || with taking into 
account U(l) gluon contributions which confirm that the 77' owes its mass to the topological 
susceptibility. 

One can say that the investigation of the internal structure of the neutral pseudoscalar 
(also vector) mesons is of great interest. As well known that very good probe of hadronic 
structure is the photon, for example, the study of electromagnetic form factors gives 
us the knowledge on the internal structure of charge particles. However, the geniune 
neutral particles like 7r°, 77, rj mesons have no the electromagnetic form factors because 
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of C-invariance. Therefore, it is needed to study the two-photon processes like P — > 77 
(P = 7T°, 77, 77') with one or both photons being off mass shell to probe their structure. 

One has to remark that theoretically the decay tt° — > 77 with both of photons being 
on mass shell is defined by the Adler-Bell-Jackiw (ABJ) triangle anomaly whereas the 
gluon contribution from Ua(1) anomaly should be taken into account under derivation of 
expression for 7/(77)77 coupling J/J. In many papers people do not use a such complication 
and just employ the parametrization of those amplitudes in terms of the decay constants 
F n , F v , F v i, and the 77 — 77' mixing angle which are determined by fitting the relevant exper- 
imental data || . A precise relation has been established between these phenomenological 
parameters and those of the chiral expansion ||. 

The more deep insight to the structure of neutral particle may be done by looking at 
the photon-dilepton decays P — ► ^l + l~ (and similar decays of the neutral vector mesons 
V — > Pl + l~ with V being p°, u, <fi) (see, review [10]). These processes are determined 



by the transition form factor FpvviQ 2 ) i n the time-like region of momentum transfer 
squared. Obviously, it is not enough to use the effective chiral lagrangians to obtain such 
form factor. It is needed to involve some assumptions on dynamics of such transitions 
like the vector dominance model (VDM) or quark models. The transition form factors 
in space-like region (Q 2 = —q 2 > 0) have been studied experimentally in the reaction 
e + e~ — > e + e~P and are planning to study in the photo and electroproduction of 77' from 
nucleons and nuclei on CEBAF's energy using the "virtual" Primakoff effect |TT[ . 



Another processes which help to shed light on the structure of the 77(77') mesons are 
decays 77(77') — > 7r + 7r _ 7 which are described theoretically by both the box-anomaly and 
resonance contribution.. The phenomenolgical analysis of these processes have been done 
in paper 
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To summarize this brief sketch, we would like to emphasize once more that many efforts 
have been devoted to describing the low-energy hadron physics in terms of effective chiral 
Lagrangians which include the Wess-Zumino anomaly |13|-||18||. There low-energy QCD 
is described in terms of mesonic degrees of freedom instead of those quarks and gluons 
with taking into account the underlying symmetries and anomalies of QCD. The hadronic 
verteces in the effective Lagrangians are structureless therefore the simple way to take into 
account the momentum dependence has been chosen by incorporating the vector mesons 



2 



either as massive Yang-Mills bosons or as the dynamical gauge bosons of a hidden local 
symmetry in the nonlinear chiral Lagrangians ||15|| , fT6|| . 

We also note that the popular methods of QCD sum rules and lattice QCD are lim- 
ited for timelike region. There are a few QCD-motivated quark models which incorpo- 
rate quark degrees of freedom and allow one to describe many-particle processes both at 
spacelike and timelike regions. One of them is based on the QCD Dyson-Schwinger 



equations with solution describing a confined quark, and another one p0| , [pT| is based on 
the phenomenological treatment of the mesons as qq states with the confinement ansatz 
to ensure the absence of quark production thresholds in the S-matrix elements of physical 
processes. Both of them have been successfully applied to describing the low-energy pro- 
cesses such as the reaction 77r — > tttt. The salient feature of our confined quark approach 



20| , [pHl is that all mesons (pseudoscalar, vector, scalar, etc.) are described as qq states by 
writing down the interaction Lagrangians and employing the compositeness condition to 
determine the meson-quark vertex coupling constants. The S-matrix elements are calcu- 
lated in the lower order in the 1/N C expansion which automatically includes the diagrams 
with intermediate resonances. 

In the present article we report on a calculation of the transition form factors of the 
processes P — > and V —>■ Pl + l~ which are related to the study of momentum 

dependence of triangle anomaly, and the f](r]') — > 77r + 7r~-decays which are related to 
the box anomaly within a confined quark model. This is a continuation of study the 



momentum dependence of anomalous processes which we started in the paper (22| . 

2 Model 

A confined quark model is specified by the interaction Lagrangian describing the transition 
mesons into quarks 

L I (x) = ig P q(x)P{x)rfq(x) + g v q(x)V^(x)^q(x) + eA il (x)q(x)Q^q(x) (1) 
with the lowest lying pseudoscalar and vector nonet meson matrices being equal to 
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The physical fields rj, rf, and lu, <ft are defined in the standard manner by introducing the 
mixing angles as 



rj = —rji sin 6p + i] 8 cos Op r{ = r\\ cos Op + r] 8 sin P 



-Ui sin Qy + ujg cos V uj = uii cos V + u 8 sin 0\ 



so that their interaction with quarks may be written as 



(2) 



—ig v T]\ -j= smop + s7 s cos dp j 

v2 

u-f 5 u + d^d _ 5 . 

+" l 9r l '' l l | 7= cos Op — S7 s sm op} 

V2 



(3) 



, .. , ir^w + d^d _ 

— n,rh^{ — sm dy + S'y^'S COSOy} 

V2 



uYu + d^d _ , 
r= cosoy — S7 P ssmoy| 



v/2 



(4) 



where S = — Oj. The "ideal" mixing angle is equal to 9j = arctan(l/v / 2) = 35.3°. 

The mixing angles are usually determined from the mass formular (quadratic or linear) 

|2gj : ^quad = ^Hn = _ 2?> o ^ ^ gquad = ggo^ Q \m = ggo Thig means ^ ^ j g yery 

close to i.e. ~ ss and we will use this value in what follows. 
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The Veneziano's model [|J which incorporates the expicit Ua{1) anomaly gives Op = 
-18°. 

In our model we use the experimental values of all hadron masses under calculations 
of matrix elements and consider the mixing angles defined by Eq. (Q) as the adjustable 
parameters. We determine 6p by fitting the available experimental data. 

The meson-quark coupling constants gu are defined by the compositeness condition 
meaning that the renormalization constants of the meson fields are equal to zero 



'M 



1 - g 2 M W M (m 2 u ) = 0. (5) 



Here U' M is the derivative of the meson mass operator and tum is the physical meson 
mass. In other words, the equation (|^) provides the right normalization of the charge 
form factor F M (0) = 1. This could be readily seen from the Ward identity 

91KAP 2 ) = 9m^^0^ = 9 2 M^P^(p,p) = F„(0) = 1. 

where Tj^{p,p') is the three-point function describing charge form factor. 
Mesonic interactions in the QCM are defined by the closed quark loops 



da v ti 



M(x 1 )S v (x 1 - x 2 )...M(x n )S v (x n - Xi 



(6) 



Here. 



is a quark propagator with the scale parameter A characterizing the size of the confine- 
ment, and the measure da v , which is essential for quark confinement, is defined to provide 
the absence of singularities in Eq. (^|) corresponding to the physical quark production: 

0(7 " G(z) = a(-z 2 )+zb(-z 2 ). (8) 



v — z 

The shapes of the confinement functions a(u) and b(u), and the scale parameter A 
have been determined from the best model description of data in low-energy processes: 

a(u) = 2 exp(-u 2 - u) b{u) = 2 exp(-w 2 + 0.4m) A = 460MeV, (9) 
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which describe various basic constants quite well (see, p0| , pl| ). 

The analytical expressions and numerical values for the meson-quark coupling con- 
stants from Eq. (||) are written down 

3.37 7T 



9P 27r V3i?pp(/4) { ^65 i 



3.24 p 



gv = 2iT\ — — = < 3.23 to 



Rvv(^ 



3.17 



We use the notation /i = m/A throughout. The structural integrals Rpp(x) and 
Rw( x ) are given in the Table I. 

Also we have the possibility to investigate the consequences of the explicit Ua{\) 
breaking in our approach. To do it, we may just assume that the octet and singlet 
coupling constants in the Lagrangian ([!]) are not equal each other 

L 81 = 4=<?(Wa 8 + £ir/iA°Vg (10) 



_ V2" 
where A = J2/3I. Using this Lagrangian in the compositeness condition ([5]) and defining 
the effective couplings of physical fields with quarks give 



2 2 2 a i 2 • 2 n 

g v = g 8 cos 9 + g x sin 6= - 



= y 8 cus u -|- y 1 sin o = — 

L P\ 

(11) 

2 2 • 2 /i i 2 2/1 1 

g r{ = g 8 sm 9 + g x cos g = 

Then we introduce two new angles characterizing the 77 — r/- mixing with taking into 
account the f/ J 4(l)-breaking: 

sin#„ = — sin# sin^„/ = — sin^. 

g v gv 

Thus the explicit £/a(1) breaking gave us two different mixing angles which are determined 
by the 77 and 7/ masses through the compositeness condition and the initial mixing angle 
9. Their numerical values are not much differ from this angle, for example, if 9 = —11° 
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then 9 V = -9.9° and 9 V > = -12.2° , if 9 = -18° then 6 V = -16° and 9^ = -20.2°. The 
available experimental data on the ?7(?7')-decays have too large uncertainties to distinguish 
these possibilities. So we will not consider the explicit C/4(l)-breaking for the time being. 

3 Radiative decays of neutral pseudoscalar and vec- 
tor mesons 

3.1 P -> 7 /+/ and V -> Pl + l decays 

The study of electromagnetic structure of neutral mesons (P = tt°, r], rf, V = p°, u, 
0) may be done in the two-photon processes: P — > 77, P — > P — > P/ _ , or in 

the processes with a change of C-parity of the mesons in the initial and final states: 
V -> PZ+Z~ and P -> W+Z". 

First, we give the necessary model- independent formulas (invariant matrix elements 
and differential distributions) which we use in our calculations. 



(a) P -> 7 Z+Z" 



1 - 



M(P^ 7 Z+r) = e 3 G P , n (m 2 P ,0,q 2 )e^ aia2 e u q^q a2 —lYl 



dT(P -> 7 /+r) 

cig 2 • T(P -»• 77) 



2a 1 



4m? 
1 ^ 



1 + 



2m 2 



l^p(g 2 )| 2 



;i2) 



Pp(<r) = ^^,0,^)7^^,0,0) 



r(p - 77) 



an 



m%Gr Prr (rrip, 0,0), 



Here p (jo 2 = mp), gi (g 2 = 0), and g are the four-momentum of pseudoscalar, photon, 
and lepton pair, respectively, and e u is the polarization vector of outgoing photon. 



(b) V -> Pl + l~. 
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m(v -> pi + r) 

dT{V -> PV~l~) 
dq 2 ■ T(V -> 7P) 37rg^ i_ ~g 



1 r 



e 2 /^ I'm 2 -m 2 /, 2 "\c-M^ a l«2 ,^„«l „«2 7_,/ti7 

e (jyp^yrriy.mp.q )e e q 1 q — J7 r J 



a 1 



4m? 



1 + 



3/2 



2m? 



(13) 



(my — mp)' 



(my + mp) ; 



3/2 



Fvp(q 2 )\ 2 



F V p(q ) = G VPl (m V} m P ,q )/G VPl (m V} m P} 0) 



a 



r(V-P 7 ) = ^ 



1 - 



ml 



mf 



Gy Pl (my,m 2 P ,0) 



Here p (p 2 = m y ), g x (g 2 = mp), and g are the four- momentum of vector, pseudoscalar 
meson, and lepton pair, respectively, e u is the polarization vector of incoming vector 
meson. 

The transition form factors F P (q 2 ) and Fy P (q 2 ) define the internal structure of the 
neutral mesons and should be calculated dynamically. 

In our approach they are defined by the diagrams in Fig. 1 for the decay P — ► , yl + l~ 
and similar diagrams for the decay V — > Pl + l~ . 

We have 



Gp 77 (p 2 ,g 2 ,g 2 ) = g Pll (p 2 ,q 2 ,q 2 ) + q 2 ^2g V p 1 (p 2 ,q 2 ,q 2 )gv 1 (q 2 )D v (q 2 ) (14) 
Gvp^(p 2 ,ql,q 2 ) = gvp^(p 2 ,ql,q 2 ) + q 2 Y,9vpv'(p 2 ,q 2 ,q 2 )gv' 1 (q 2 )Dy(q 2 ) (15) 

V 

Here the first term corresponds to the triangle diagram whereas the second to the vector 
meson pole diagram in Fig. 1. One should emphasize that in our approach the con- 
tribution from the pole diagram vanishes if the photon is on mass shell (q 2 = 0) as a 
consequence of gauge invariance. 



The quantities appearing Eqs. (14-15) are written down 
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9p^p 2 ,qlql) = ^- 2 jC Pll R PVV (p 2 /A 2 ,ql/A 2 ,ql/A 2 ) 

We do not write the expressions for gvpv because of we will consider the only decay 
uj -> In this case g^o p = (3g w g p g OJ /2V2n 2 )R PVV (p 2 /A 2 ,qf/A 2 ,q^/A 2 ). 

The structural integrals and SU(3)-factors entering these expressions are given in Table 
I and II. 

The Eqs. ( [L4} - J15D may be written in the forms which are more convenient for compar- 
ison with the VDM approach: 



G 7T o^(p 2 ,qf,q 2 ) 



g^(p 2 ,q 2 ,q 2 )^l + q 2 



\cM 2 ) + \cM 2 ) 



o sin cos5p 



3 sin 5 P 



G vljy (p 2 ,qf,q 2 ) = g v >y~ / (p 2 ,q 2 ,q 2 )-{l + q 2 



5 cos <5p — v2 sin <5p 



where 



3cos5 P 2 cos5 P 2 v^sin^p 2 
2 C M ) + — g~ ) 3 C^q 



(16) 



with being a vector meson propagator. 

The expression for the renormalized vector propagator obtained within our model [24 
in one-loop approximation is written 
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m 2 v R v {m 2 v /h 2 )-q 2 R v {q 2 /h 2 ) 
where the obvious equality (see, Table 1) (xRy(x))' = Ryy(x) provides the residue at 
q 2 = my to be equal to one. 

The imaginary part of the p(o>)-meson propagators is not very important for 7r and 
77 decays but it should be taken into account somehow for rf decay. In our approach 
the imaginary part should appear as a result of summing up the two- loop diagrams. It 
will be proportional to the whole decay width of a vector meson, for example, for p- 
meson it reads T p {q 2 ) = m p {g 2 n7T {q 2 )/A87i)[l - 4m 2 n / 'q 2 ]^ 2 Q{q 2 - Ami) wh ere g p7T7T {q 2 ) = 
(3glg p /4V2n 2 )R VPP (q 2 ). 

Then the vector propagator may be written as 

D („2x Ryy(m 2 v /A 2 ) 

V[q ) m 2 v R v (m 2 v /A 2 ) - q 2 R v {q 2 /h 2 ) - im v T v R vv {m 2 v ) 1 ] 

with Yy being the experimental value for vector meson decay. We drop out the q 2 de- 
pendence of decay width since it is the resonance region. The extra factor Ryvijny) in 
denominator is introduced to provide the correct normalization of propagator at q 2 = my. 
Also we will use the same magnitude of width for p and 00 mesons T v = 151 MeV because 
of uj has too small width (about 8 MeV) to be seen in the experiments and its account 
theoretically will just give unpleasant extra pick. Of course, one drops 0-width because 
it is not important at energy less 1 GeV. 

For comparison, we will use the simplified VDM-form for vector propagators 

= <-/-*MV (19) 

We will discuss the numerical results for decay widths and the form factors in the final 
subsection. 

3.2 77(77') — > 7r + 7r~7 decay- 
First, we give the model-independent expessions for invariant matrix element decay width 
and distributions of the P — > 7r + 7T~7-decay (P=?7 or 77'): 
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M"(P - 7T+7T-7) = ee^plp^Gp^-^S!, s 2 , s 3 ) 

(20) 



r(P^7T + 7T~7) = ^3" / dS2S<2 j ds l( S l - S l)( S l- Si)\G P ^ n + n - ry (s 1 ,S 2: S 3 ) 



Here, the kinematical variables Sj are defined as 

si = (g + p+) 2 s 2 = (p_+p+) 2 s 3 = (g + p_) 2 , si + s 2 + s 3 = m P + 2m 2 . 
so that they vary in the intervals 4m 2 r < s 2 < mp and sj" < si < s^" where 



sf = + (1/2) (m 2 , - S 2 )[l ± ^1 - 4mJ/s 2 ]. 

In our approach the matrix element of the decay P — > 7r + 7r~7 is defined by the 
diagrams in Fig. 2. The expression for G P ^ 7T + 7T - 1 is written down 



Cp_ >7r + 7r - 7 (si, s 2 , S3) — <7p (si, s 2 , s 3 ) + 2#p p7 (s 2 )# p7r7r (s 2 ).D p (s 2 ) 

g Pfn (s 2 ) = C P j^9pR PVV (m 2 P /A\s 2 /A 2 ,0) 

9^-^,8,, s 3 ) = C P j- 3 ^^R n ( Sl /A 2 ,s 2 /A 2 ,s 3 /A 2 ) 

R n ( Sl ,s 2 ,s 3 ) = f '^ a 5{l-Y t a i ){-a'i < -Dl^) + a\-Dp-)-a\-Dr + )} 
J i=i 

where 

D± +1 = a^a^mp + a 2 (ai + o 3 )m 2 + a 1 a 3 s 2 + o 2 o 4 Si 
D 4 7 " = a 2 a 3 m 2 P + (0:30:4 + oio 2 )m 2 + 01O3S3 + o 2 o 4 si 
Dj~ + = aia 2 m P + o 3 (o 2 + o 4 )m 2 + Oi0 3 s 3 + o 2 o 4 s 2 

The group coefficient C P is equal to — sin^p for r\ and cos^p for r]'. 

One trivially finds that Pn(0, 0, 0) = 1/3. In the physical region the box-contribution 
varies very slowly around this value so that we may approximate R n by this value with 
quite good accuracy in the calculations. 
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This gives the following formula for the decay width 



mp 

T(P - vr+vr-7) = ^ J dmE*I#\G®l^^{m*)\* . (21) 

2m, 



with E n = {mp — m 2 )/2mp and P T = (l/2)ym 2 — Am 2 being the photon energy in the 
ingoing meson rest system and the pion momentum in the dipion rest system, respectively. 
The label (0) upstairs G means that we use Rq = 1/3 in the whole physical region. 
The effective mass spectrum of the 7r + 7r~ system is defined as 

' /r " E*P*\A(™l«)\ 2 - (22) 



dm™ 12tt 2 7 *' ™ 
The photon energy spectrum is defined as 



" ' E^\G^ + (s 2 )\ 2 . (23) 



with S2 = mp{mp — 2E 1 ). 

We will discuss the numerical results in the final subsection. 

3.3 Discussion of the numerical results 

The dependence of the ratios of the theoretical values of decay widths to their experimental 
averages on the mixing angle are shown in Fig. 3a Fig. 3b for the vector propagator shape 
described by Eq. (0) and Eq. (|T9|) , respectively. Visually, one can see that the value of 
the mixing angle 6p = —16.5° in Fig. 3a seems may give quite resonable coincidence with 
the experimental data except the decay i] — > 7r7r7 which theoretical value is 1.5 times 
larger than the average experimental one. One has to note that the value 9p = —16.5° is 
supported by two recent theoretical analysis: (1) r],r]' — > 77* and D s — > Vjlvjv(lv which 
gives 6p = —16.7° ± 2.8° p5| , and (2) J/ifj decays into a vector and pseudoscalar meson 



which gives 9 P = -16.9° ± 1.7° 

First, we discuss the decays T](j]') — > 7/i + /i _ and uo — > 7r° / u + / u _ . The transition form 
factors decsribed the observable spectra Eqs. (P^TB|) have been parametrized in the pole 
approximation |TUj as 
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F ^ 2 ) = _ T7U A, = 0.72 ±0.09 GeV, (24) 

Q / r; 

F^{q 2 ) = ; Vttt A w = 0.65 ± 0.03 GeV. (25) 

1 - q 2 /Ai 



As it is seen from Eqs. (J24]-|25D a value of the characteristic mass, A, agrees well with 
vector dominance (m p = 0.77 GeV) for //-decay, and differs from those by four standard 
deviations for cu-decay. Our form factor (see, Fig. 4) goes slightly lower than the curve 
fitted the experimental data but it is still within the experimental uncertainties. The 
behavior of u — 7r-form factor (see, Fig. 5) almost coincides with the fitting curve. 

The behavior of the form factor describing the rf — ► 7/i + p _ decay is more interesting 
since the p° and u poles occur in the physical region allowed for the muon pair spectrum. 
Experimental data on the r/-meson transition form factor, the VDM predictions with 
taking into account the p-meson decay width in the vector propagators, and the results 
of our calculations are shown in Fig. 6. A more detailed analysis is difficult to perform 
because of poor statistics. 

Now let us discuss the decay rf — > 77r + 7r~ and related process rf — > p°7. There are lots 



of experiments devoted to the study of these decays ||29||-|39|. Most of them are looking 



at the effective mass spectrum for the 7i + n system in the decay rf — > 7i + n 7 (see, Fig. 



7). It was shown in [39| that an attempt to fit the distribution in the framework of the 
cascade decay model rf — > p°7, p° — ► 7r + 7r~ with simple choice of the vector propagator 



dN ElPlm* 

'J 7T 7T7T 

OC 



dm^ (m^ - m 2 p ) 2 + m 2 p T 2 p 
with T p being the "dynamical" width of p-meson is not successful because of the distri- 
bution obtained in experiment is much harder than the theoretical one. It was needed 
to assume that either the propagator form should be modified significantly without res- 
onable theoretical justifications or a nonresonance process gives significant contribution 
to this decay. The last assumption was realized |39| by modelling the distribution by the 
following form 

dN 1^ ^^ + fe<a " 



dm n7T ~ 7 * 7r7r ml n -m 2 p + im p T p m\ 
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As a result of fitting the experimental spectrum it was found that £ = 2.78 ± 0.46 and 
a = -1.07 ±0.08. 

Recent measurements of the it + tt~ mass spectrum in the decay rf — > 7r + 7r _ 7 with 
the Crystal Barrel detector also confirmed the existence of a non-resonant coupling 
7r + 7r~7. We plot their data on Fig. 7a and compare them with our results and VDM 
prediction. It is readily seen that our results coincide quite well with the experimental 
data whereas the VDM predictions are significantly above them up to the resonant region. 

The existence of the non-resonant contribution clearly comes from the existence of the 
box anomaly in chiral theory. One of the attempts to incorporate vector mesons in the 



chiral Lagrangian has been done in |16| . It was found that the amplitude of the decay 
7] — > 7r + 7r~7 is given by 



M(r] 



7) 



An 2 F2 



2 m 2 - p 2 



V3£ 



cos 9 



I21_ 
3/o 



sin 9} (26) 



where /§ = 1.25F n and /o = 1.04i^. A numerical result for the decay width obtained with 
9 = —20.6 was T(r] — > 7r + 7r~7) = 62 eV which is close to the experiment r expt = 58 ± 6 
eV. The direct generalization of this model for r/-decay by adding the imaginary part like 
irripTp to the denominator of propagator and replacing cos# to sin# and sin^ to -cos 9, 
respectively, gives the spectrum and the decay width listed in Fig. 7b in a agreement 
with the experimental data and our results. 

Now we discuss the rj — > 7r + 7r~7-decay. From the theoretical point of view this is 
simpler than the 7/-decay since it is up to the resonance region. But experimental data 
are very poor, actually, there are only two experimental papers discussing this issue PS] 
and [2~7| where the 7-ray energy spectrum was measured. To compare the theoretical 
results with the experimental data, two forms of function f(s) entering to the matrix 
element M(r] — > ti + ti^j) oc / \s)e^ ual3 e^p^p+p^ have been used |28|: f(s) = 1 and f(s) = 
m 2 /(m 2 — s). The efficiency used in this procedure has not been given in explicit form. 



We are able to reproduce the results of the paper [28| with the efficiency function chosen 
in the form 



eff(£ 7 ) = EZ{E&> - E 7 ) 



(27) 
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where a — b — 1.1. Using this efficiency function we recalculate the experimental data 
according to 

corrected data = data/eff (£" 7 ) 

and put them in Fig. 7a and Fig. 7b. One can see that the theoretical results from 
QCM, VDM, and the chiral approach implemented VDM [[16] coincide each other and fit 
very well the experimental data. Finally we summarize our predictions for the rates of 
radiative decays of r] and rj' mesons in Table III. 
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Table I. Variuos functions entering the QCM calculations of meson vertecies. Functions 
a(u) and b(u) characterize the confined quark propagator Eq. (|p-|9|). 



Rpp{x) 


= #o + 




txyy \X ) 


= #o + 




Kv{X) = 


A) + f 


Jdub(-f)(l + u/2)y/l-u 


R V pp{x / 


= £o + ! Jdub{-f)Vl-u 


Rpvv{%- 


y,z) = 


3 

/ d s a5(l — Y<oii)a{—aia2X — a\a^y — a^a^) 


Rpvv{%- 


y,o) = 


(xR PVV (x, 0, 0) - yR PVV (y, 0, 0))/(x - y) 


Rpvv{x. 


0,0) = 


\Sdua( T)in(;:^) 


CO 

Aq = J dua(u) = 


OO 

= 1.09 B Q = J dub{u) = 2.26 
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Table II. The SU (3 )-f actors entering the calculations of P — > 77 (P — > 7/+/ ); 
V -> P7 (V -> PZ + Z~), P -> F 7 ; and V -> 7 -> Z+Z~ decays. 



C P77 = 3tr(A p Q 2 ) 


P 


1 




-(5 sin<5p + \/2 cos<5p)/3 


V 


(5cos5 P - y/2waS P )/3 


7/ 



C VP , = tr(Q{\ p ,\ v }) 


\/P 7 


2/3 


P7T 7 


2 


dJ7T 7 


—2 sin5p 




-(2/3) sin 8 P 


CJ777 


-(4/3) cos 5 P 




2 cos Sp 


pn'i 


(2/3) cos 5p 




-(4/3) sin 5p 


07?'7 



Cy 7 = tr(A y Q) 


V 


1 


P° 


1/3 




v^/3 
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Table III. Comparison of our predictions for radiative decay widths of neutral 
pseudoscalar and vector mesons with the available experimental data. Results are given 

for the 1] — r/-mixing angle Op = —16.5°. 





QCM 


Expt. 


7T° — > 77, eV 


7.2 


7.7 ±0.6 


r\ — > 77, KeV 


0.45 


0.46 ±0.04 


v' ~ ¥ 77? KeV 


4.2 


4.26 ±0.19 


p° — > 777, KeV 


38 


57 ± 12 


77 — > // + //~7, eV 


0.35 


0.37 ±0.06 


77 — > 7r + 7r _ 7, eV 


92 


58 ±6 


77' — > // + //~7, KeV 


0.014 


0.021 ±0.007 


7/ — > 7r + 7r _ 7, KeV 


47 


56 ±9 


77' -> p° 7 , KeV 


38 


61 ±8 


a; — > 7T7, KeV 


590 


716 ±51 


a; — > tth + h~, KeV 


0.49 


0.81 ±21 
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Fig. 1. The P — > decay. 
Fig. 2. The P — > 7r + 7r _ 7 decay. 

Fig. 3. The dependence of ratios of theoretical values for the radiative decays of r\ and rj' 

mesons to their experimental averages on the mixing angles calculated with the shape of 

vector meson propagators defined by (a) Eq. fllffi), and (b) Eq. (|Hf). 

Fig. 4. The form factor of rj — > / jfi + fi~ decay. 

Fig. 5. The form factor of uo — > 7r° / u + / u _ decay. 

Fig. 6. The form factor of v( — > 7yU + /i~ decay. 

Fig. 7. The mass spectrum of rj' 77r + 7r~ decay: 



(a) QCM and VDM (b) Chiral model |T§ and VDM. 
Fig. 8. The photon energy spectrum of rj — > 77r + 7r _ decay: 



(a) QCM and VDM, (b) Chiral model 0] and VDM. 
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